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Abstract 

We investigate — as an alternative to usual Monte Carlo Renormalization 
Group methods — the feasibility of extracting nonperturbative QCD j3 func- 
tions directly from a lattice analysis of correlations between the action and 
Wilson loops. We test this correlation technique numerically in four dimen- 
sional SU(2) gauge theory, on a 16 4 lattice at f3 = 2.5 and find very promising 
results. 
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I. INTRODUCTION 



Determinations of QCD /3 functions beyond perturbation theory are becoming more and 
more important. This holds in particular for the lattice laboratory, where the quality of 
data allows for increasingly detailed physics investigations, like scaling analyses and ther- 
modynamic studies in anisotropic lattices. 

Much work has been done in the past to extract (5 functions from the response of lattice 
observables to changes in the parameters (coupling and, in the case of full QCD, quark 
masses) [1-4]. In these analyses, the j3 function is typically retrieved from finite difference 
measurements, at fixed values of physical observables (like physical hadron masses). In 
practice one has to operate at sufficiently different values of the parameters, however, i.e. 
not really locally. 

It is well known in statistical mechanics that response functions are related to correlations 
and therefore can be determined without recourse to such finite difference methods. In 
lattice gauge theories this correspondence has become known in the context of lattice sum 
rules [5-7]. The latter have been considered in the past as mere consistency checks for the 
underlying (action and energy) distribution functions, assuming the j3 function to be known. 
Given the recent progress in lattice technology, however, we think it is timely to turn things 
around and exploit sum rules for a nonperturbative determination of the j3 -function. 

In this paper, we will test this approach. For simplicity we will work with SU(2) gauge 
theory. We shall extract the (5 function, Bp = 8(3/ d In a, from the action-density in presence 
of a static quark-antiquark pair, ie. from correlators between (smeared) Wilson loops and 
plaquettes. The method can readily be generalized to other response functions, such as 
B K = dn/dlna, appearing in full QCD. 

In the next section we will introduce the idea of the method for the case of dynamical 
fermions. This includes of course the pure gauge theory case that we pursue in the subsequent 
numerical part. The details of the numerical techniques are explained in section 3. The 
results are compared to results from interpolation methods in the concluding section 4. 
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II. THE GIST OF THE METHOD 



We consider the (isotropic) lattice action 

S = S g (/3;U) + S f (K;U,^) , (1) 

where S g (Sf) denotes the gluonic (fermionic) part. The most common form of S g is the 
Wilson action 

S g = (3S W = (3 (2) 

with 

S, v {n) = 1 - ^Re Tr (U,(n)U v (n + /})[/> + v)Ut(n)) (3) 

where N denotes the number of colors. 
The fermionic action is given by: 

Sf = ^i) m M mn il) n (4) 

m,n 

with the fermionic matrix 

M = 1 - kD. (5) 

tp and tp are to be understood as vectors in color, Dirac, and position spaces, M as a 
ANV xANV matrix. The hopping term, D, corresponds to the tpd^tp part of the continuum 
Dirac action. The results presented here are quite general and do not depend on the specific 
form of the action. 

Let us now illustrate the basic idea for computing the j3 functions. We consider the 
expectation value of a (gluonic) operator, O, 

^ ' ~ J dUdiPdiPe- s ( u ^ ' ^' 
and compute the responses under changes in the couplings j5 and n: 
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r\o]- dln{0) - ( {0Sa) IS)) (7) 

and 

^ = ^? = W-<w) • 

Note that the chiral condensate, (tptp) = (TrM -1 ), can be obtained from inversion of the 
fermionic matrix on fixed gluonic background fields. 

The expectation values of S g and tptp are to be understood without normalization, i.e. 
without division by factors 6V or 4NV. Eq. (8) has been obtained from the identity 

(^-) = -n{ijDij) = (tpMtp) - = ^NV - (ijip) . (9) 
\ a In k J 

If the derivatives of InO in respect to the couplings can be related to derivatives of 
dimensionful quantities, the right hand sides (RHS's) of Eqs. (7) and (8) can be measured 
to determine the latter. Let us consider as our test operator O a Wilson loop 
encircling a (rectangular) area A. In the confined phase one obtains the area law 

W(A) oc exp(-KA) , (10) 

where A has to be kept sufficiently large to suppress perimeter terms and (in case of dynami- 
cal fermions) sufficiently small for string breaking to be ineffective. The constant K is related 
to the physical string tension, a = K(/3, «)a _2 (/5, k). This relation allows to determine the 
lattice resolution a. 

From Eq. (10), one obtains: 

ld\nW(A) da 2 dlna , x 

A a In 7 a In 7 am 7 

where 7 = (5 or k. Finally we find from Eqs. (7) and (8), 
and 
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B K « -2AKkC- x \W{A)\. (13) 

Note, that the relations Eqs. (11) and (12) have been derived by Karsch long ago [8], in the 
context of thermodynamics. 



III. NUMERICAL ANALYSIS 

Though, in theory, Eqs. (12) and (13) offer the possibility for a direct measurement of 
the j5 functions, in practice, strong fluctuations of the correlators, in particular for large 
Wilson loops, might affect the applicability of this method. We can demonstrate, however, 
that appropriate use of smearing and noise reduction techniques in conjunction with the ex- 
ploitation of the known subasymptotic form of Wilson loops (from transfer matrix methods) 
leads to substantial improvement in correlation signals, and provides success. 

We start from the spectral decomposition of a rectangular Wilson loop with a spatial 
extent of R and a temporal extent of T lattice units: 

W(R, T) = C(R)e~ v ^ T (l + Y, f|§ e ~ A W ) > (14) 

where AVi(R) = Vi(R) — V(R) are gaps between ground state and excited state potentials. 
The overlaps C(R) and Ci(R) are positive, with their sum normalized to one (C(R) + 
J2i Ci(R) = 1). Manipulation of the spatial parts of the Wilson loop leaves the eigenvalues 
V and Vi unaffected but alters the amplitudes C and C;. 

Let us differentiate a smeared Wilson loop (with fixed lattice R, T) with respect to j3: 

^ m R, T)] ^J^^ (15) 

_ d\nC(R) _ ( dV(R) d(R) dAV^R) r -AvuR)r\ , h fi . 

~ d(3 { d(3 + C(R) d(3 ) + ' ' ' ' 1 ' 

Ground state enhancement implies an additional suppression of the (exponential) correction 
terms to the linear behaviour, which is welcome. Therefore, we construe a combination of 
suitable spatial paths to deplete excited states contributions up to a few percent, by use of 
APE-like smearing methods [9,10]. 
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As a side product to our previous work on flux distributions [10], we have carried out 
simulations on a 16 4 lattice at j3 = 2.5, with 8680 configurations 1 . In Fig. 1, we illustrate 
the numerical situation for R = 3. We observe a clean linear behaviour of Cp(T), starting 
out from T as low as T min = 2, which is indicative for the suppression of excited states. 
Note that the slope is - apart from the normalization 2K - our signal for — dV a ^ ■ 

After having extracted dV(R)/d(3 from linear fits in T > T m ; n , starting from reasonable 
values of T min , we separate the potential into a self energy part from the static propagator 
— which will diverge in the continuum limit — and a physical part, 

V(R) = V + V ph (R) = V + av(aR). (17) 

v(r) is the potential, measured in physical units, and does not depend on the lattice spacing. 
Starting from this separation ansatz, one finds 2 

By applying a two parameter fit to the RHS in accord with the expected behavior, Bp 
can be extracted. 

In the confined phase of QCD, the potential can be fitted to the ansatz V p h = KR — e/R 
with high accuracy, starting from Ra 0.3 fm. From this parametrization, one expects the 
behaviour 

with a constant c that amounts to 1/(2K) times the derivative of Vq with respect to (5. It 
is nice to see at this point, that the Coulomb like term has disappeared, due to its scale 



For the present feasibility study we have not attempted to average over all possible orientations 
and locations of Wilson loops nor have we made use of link integration techniques. The cost of the 
method will readily be reduced by a factor 30-50, once this is done. 

2 The RV' term has been missing in the original form of the sum rules [5], as has been noticed in 
ref. [11]. 
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invariance 3 . Indeed, this will enable us to perform the analysis down to very small values of 
R. 

The actual numerical situation is depicted in Fig. 2, where the values for — ^ 9 V Q ^ are 
plotted versus R. The data indeed exhibits the linear behaviour in R predicted in Eq. (19), 
with a slope 

Bpl 2 5 = -4.03 ± 0.23 , (20) 

for -R m ; n = 3. Table I collects the results from various fits, with different cuts. We find 
the value of the (5 function to be fairly stable under reasonable changes of -R m ; n and T m ; n . 
Note that the nonperturbative result differs appreciably from the perturbative two loop 
expectation for the (5 function which reads B 1 = —2.53 for SU(2) gauge theory at (5 = 2.5. 

IV. DISCUSSION 

It is interesting to compare our local value for B 1 with previous nonperturbative results 
that are based on interpolations to a(f3). A fit to the string tension values across the region 
2.3 < p < 2.85 renders the estimate [10] B 1 m —2.94, whereas the analysis of the critical 
temperature data presented by Engels et al. [12]led to the the value B~ l m —3.21 for (5 = 2.5. 
The difference between these "global" estimates might be interpreted as a 10 % interpolation 
uncertainty. 

On the other hand we find the result from our present local method, Eq.(20) to differ 
by more than two standard deviations from the "global" estimates. A high statistics scan 
through the /5-region with the local method is obviously needed to study the situation in 
more detail and settle the issue. 

It looks rather promising to pursue this method in the case of full QCD, as knowledge 
about the K-response function will be very helpful in guiding full QCD simulations through 



3 We thank H. Dosch for discussions on this point. 
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coupling parameter space. The same holds for the finite temperature field theory with 
anisotropic couplings. 

We would like to stress that the local method will be particularly efficient in a high 
dimensional coupling space, compared to finite difference methods that combine different 
Monte Carlo runs. 

For the purpose of our feasibility study we have used the string tension, being the most 
precise dimensionful quantity accessible in lattice simulations [13]. It would be interesting 
to benchmark the power of correlation methods on other observables. 
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FIGURES 
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FIG. 1. The correlator Cp[W(R,T)]/(2(3K) versus T at (3 = 2.5 and R = 3. The slope at large 
T is proportional to the derivative of the potential in respect to (3. 
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FIG. 2. The combination — ^ versus R at f3 = 2.5. The slope at large i? is the inverse 

of the (3 function, Bp 1 . 
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TABLES 

TABLE I. Results on the (3 function B and the constant c (Eq. (19)) at (3 = 2.5 for different 
cuts T min and R min . 



T m j n 
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